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Abstract 

This paper describes a new algorithm for determining all discrete contact symmetries 
of any differential equation whose Lie contact symmetries are known. The method 
is constructive and is easy to use. It is based upon the observation that the adjoint 
action of any contact symmetry is an automorphism of the Lie algebra of generators 
of Lie contact symmetries. Consequently, all contact symmetries satisfy various com- 
patibility conditions. These conditions enable the discrete symmetries to be found 
systematically, with little effort. 

1. Introduction 

Discrete symmetries of differential equations are used in various ways. They map solutions 
to (possibly new) solutions. They may be used to create efficient numerical methods for 
the computation of solutions to boundary- value problems. Indeed, there is currently much 
research into techniques for constructing numerical methods that preserve various types of 
symmetry [2, 4, 10]. Discrete and continuous groups of symmetries determine the nature 
of bifurcations in nonlinear dynamical systems. Equivariant bifurcation theory describes 
the effects of symmetries, but it may yield misleading results unless all symmetries of the 
dynamical system are known [3, 6]. 

In general, it is straightforward to find all one-parameter Lie groups of symmetries 
of a given system, using techniques developed by Sophus Lie more than a century ago 
[1, 11, 12, 14]. Yet, until recently, no simple method for finding all discrete symmetries 
was known. Ansatz-based methods can be used to find discrete symmetries belonging to 
particular classes, e.g. [5] , but such methods cannot guarantee that all discrete symmetries 
have been found. The main difficulty is that, commonly, the determining equations for 
discrete symmetries form a highly-coupled nonlinear system. Reid and co-workers have 
developed a computer algebra package aimed at reducing this system to a differential 
Grobner basis [13], but the method is computationally intensive and seems not to have 
been widely used. 
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A new approach to the problem of finding discrete point symmetries has recently been 
described by the author [7]. Instead of trying to solve the symmetry condition directly, 
one first examines the adjoint action of an arbitrary discrete point symmetry upon the 
Lie algebra of Lie point symmetry generators. This yields a set of necessary conditions 
which simplify the problem of constructing all discrete point symmetries. Of course, 
one must know the Lie algebra of Lie point symmetry generators, which must be non- 
trivial. However, this is not a severe limitation, as most differential equations of physical 
importance have a non-trivial Lie algebra. 

Most techniques of Lie symmetry analysis readily generalize to contact symmetries, 
and the aim of the current paper is to show how all discrete contact symmetries of a given 
differential equation may be found systematically, using an extension of the algorithm for 
determining discrete point symmetries. Perhaps surprisingly, many differential equations 
without non-point Lie contact symmetries have non-point discrete contact symmetries. It 
is useful to be able to derive these (non-obvious) symmetries systematically because, like 
all symmetries, they constrain the behaviour of solutions of the differential equation. 

The new method is described first in the context of ordinary differential equations 
(ODEs) of order n > 3, then adapted to treat a finite subgroup of the infinite group of 
contact symmetries of a given second-order ODE. Finally, it is shown that the method can 
be extended to partial differential equations (PDEs). 



2. The algorithm applied to ODEs 

A diffeomorphism 

r : (x,y) i ^ (x,y) 
is a symmetry of the ODE 

y^=tu(x,y,y',...J n -V) (2.1) 

if it maps the set of solutions to itself, i.e. if 

^ (n) =w^,y',.-i M ) when (2.1) holds. (2.2) 

Here the functions y^ are obtained by prolonging the diffeomorphism T to derivatives, 
using 

*" >s t (* (0) = *)' < 2 ' 3 > 

where 

D = d x + y'd y + y"d y , + ■■■. 

For point symmetries, x and y are functions of x and y. Contact symmetries are more 
general than point symmetries, because x, y and y' are functions of x, y and y'. One- 
parameter Lie groups of contact symmetries are obtained in a similar way to Lie point 
symmetries, by linearizing the symmetry condition (2.2). Specifically, 

x = x + e£(x,y,y') + 0(e 2 ) , 

y = y + erj(x,y,y') + O (e 2 ) , 

y{k) = y(k) + erj (k) y> y', • • • , y (*)) + Q ( £ 2 ) , k > 1. 
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The contact conditions require that the O(e) terms should be expressible in terms of the 
characteristic function Q(x,y,y') = r] — y'£, as follows 

t = -Q y >, r] = Q-y'Q y ,, r/ fc > = D k Q - y {k+1) Q y ., k > 1. (2.4) 

In particular, 

v {1) ( x ,y,y') = Qx + y'Qy 

Lie point symmetries of ODEs of order n > 2 are Lie contact symmetries whose charac- 
teristic function is linear in y'. 

The set of all infinitesimal generators of Lie contact symmetries of a given ODE of 
order n > 3 forms a finite-dimensional Lie algebra, C, which can generally be determined 
systematically [14]. Given such an ODE, suppose that C has a basis 

Xi = (,i(x,y,y > )d x + 7] i (x,y,y')d y + 7] < i 1) (x,y,y')d yl , i = l,...,N, (2.5) 

where rji, r^ 1 ^ are obtained from the characteristic function Qi(x,y,y') using (2.4), 
and N = dim(£). The structure constants, c k j, for the basis (2.5) are determined by 

[X i ,X j ]=($ j X k . (2.6) 

(Summation is implied when an index occurs twice, one raised and once lowered.) The one- 
parameter Lie group of contact symmetries corresponding to a particular Xi is obtained 
by exponentiation. We use the notation 

Fi(e) : {x,y,y') ^ (e eX *x, e eX *y, e eXi y') . 

Suppose that 

r : (x, y, y') ^ ( x (x, y, y'), y(x, y, y'), y'(x, y, y')) (2.7) 

is a contact symmetry of the given ODE. Then the contact transformation obtained by 
the adjoint action of V upon Tj(e), 

f i (e) = rr i ( e )r- 1 , (2.8) 

is also a contact symmetry, for each e in some neighbourhood of zero. Therefore, for each i, 
there is a (local) one-parameter Lie group of contact symmetries 

f i(e) : (a, y, y') h-> [e^x, e^y, e eXi y'^ , (2.9) 
whose infinitesimal generator is 

X = TXiT- 1 = £i(x, y, + th(x, y, y')d § + r?f \x, y, y')dy,. (2.10) 
Consequently 

XgC, i = l,...,N. 

The generators {Xi}f =l are simply the basis generators {Xi}f =l with [x,y,y') replaced 
by (x,y,y r ). Therefore the set {Xj}^ is a basis for £, and so each Xi can be written as 
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a linear combination of the X^s. Also, the mapping i— > is an automorphism of £ 
which preserves all structure constants, i.e. 

[Xi,Xj] = c\jX k when (2.6) holds. (2.11) 

These results generalize to partial differential equations, and are summarized as follows. 

Lemma 1. Every contact symmetry T of an ordinary differential equation of order n > 3 
induces an automorphism of the Lie algebra, C, of generators of one-parameter local Lie 
groups of contact symmetries of the differential equation. For each such T, there exists a 
constant non-singular matrix (b\) such that 



Xi = b\X h (2.12) 
This automorphism preserves all structure constants. 

Lemma 1 yields the following PDEs for the unknown functions x(x, y, y') and y(x, y, y'): 

X i x = b l i X l x = b l &(x,y,y'), i = l,...,N, (2.13) 

X i y = b\X l y = b l i r }l {x,y,y'), i = l,...,N. (2.14) 

This set of 2N first-order PDEs, together with the contact condition 
dy 

y'(x,y,y') = 4jjL, 
dx 

provides necessary, but not sufficient, conditions for T to be a contact symmetry. The 
contact condition yields the following pair of PDEs, because y' is independent of y": 

Vx + y'y y = (xx + y'x y )y', y y > = Xyiy' . (2.15) 

N.B. The lemma gives a further N PDEs 

Xtf = b[X iy ' = b\r,f ] (x, y,y'), i = l,...,N, 

but these add nothing new, for they are a consequence of (2.13), (2.14) and the contact 
condition. 

To find all discrete contact symmetries, proceed as follows. First solve the system of 
PDEs (2.13), (2.14), to obtain (x,y) in terms of x,y,y',b\ and some unknown constants 
(or functions) of integration. If N is sufficiently large, it may be possible to solve this 
system algebraically (by eliminating the derivative terms); otherwise, the method of char- 
acteristics should be used. Incorporate the contact condition (2.15); this generally reduces 
the number of candidate solutions (x,y). Finally, use the symmetry condition (2.2) to 
determine which of these solutions are symmetries. The continuous symmetries may be 
factored out at a convenient point in the calculation. The remaining discrete symmetries 
are inequivalent under any continuous symmetry, and form a discrete (but not necessarily 
finite) group. 

If C is non-abelian, some of the structure constants are non-zero, enabling the matrix 
B = (pi) to be simplified before any of the above calculations are done. Substituting (2.12) 
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into (2.6), and taking (2.11) into account, we obtain the following constraints on the 
components of B: 

cYJ>W = 4 b k- (2-16) 
It is sufficient to restrict attention to equations (2.16) with i < j, because the structure 
constants are antisymmetric in the two lower indices. Moreover, at least some of the 
continuous symmetries can be factored out using their adjoint action upon the generators 
in C, (see [11]), 

Ad(exp(e^))^ = X t - e,LY„A^ + ±[Xj, [X„X t ]] -■■■ = a\ ( £j , j)X p . (2.17) 

Let A(j) denote the matrix whose components are a?(ej, j), as defined by (2.17). The 
system (2.12) is equivalent, under the group generated by Xj, to 

Xi = b\Xi, 

where b\ are the components of 

B = A(j)B. (2.18) 

The mapping B i-> B does not affect (2.16), so we will drop tildes as soon as each 
equivalence transformation has been made. Each generator Xj is used in turn to simplify 
the form of B. 

To illustrate this procedure, consider the two-dimensional non-abelian Lie algebra o(l), 
with a basis {Xi, X2} such that 

[Xi,X 2 ]=Xi. (2.19) 

The only non-zero structure constants are 

1. 



c 12 



c 21 



Therefore (2.16) gives 
b\bl-b\b\ = bl 
and hence 

0' 



B 



bl 



1 



The matrices representing the adjoint action of the continuous group on the Lie algebra 
are 



A(l) = 



1 

-ei 1 



A(2) = 



e ea 
1 



b 1 

Applying (2.18), first with j = 1, e\ = -jj, then with j = 2, e 2 = — In \b\\, we obtain 



B 



a 
1 



where a G {—1, 1}. 



(2.20) 



The reduced form of the matrix B is specific to this particular Lie algebra, and is inde- 
pendent of the ODE whose Lie point symmetries are generated by the algebra. 

Simplified matrices for other non-abelian Lie algebras can be found by the same tech- 
nique. However, if C is abelian, the entries of B cannot be determined a priori. 
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3. Examples 

The third-order ODE 



v" 2 v" 

y'" = y —- v - 

x y' 



(3.1) 



has a two-dimensional Lie algebra of generators of Lie contact symmetries. These Lie 
symmetries are actually point symmetries, and C is isomorphic to a(l). The basis of C, 

X 1 = d y , X 2 = |<9 X + yd y + V ^d y/ , 

has the commutation relations (2.19), and therefore B is given by (2.20). The system of 
PDEs (2.13), (2.14) amounts to 



~X lX 


Xif 




'a 0" 


_X 2 x 


X 2 y _ 




1 



" 


1 " 




' 


a 


X 






X 




- 2 


y 




- 2 


y 



The general solution of this system is 

x = xp(t), y = ay + x 2 q(t), where t = — . 
The contact condition gives 
y' = xr(t), 

where 

at + 2q — tq = (p — tp)r and q = pr. 

(Here a dot over a function denotes its derivative with respect to t.) Re-arranging these 
conditions, we obtain 



and the contact condition is satisfied if and only if 
pr — pr = a. 



(3.2) 
(3.3) 



It is convenient to work in terms of p and r, because the ODE is invariant under translations 
in y, and so y does not occur in the symmetry condition. The prolongation to second and 
third derivatives is 



y 



y 



r + (y" — t)r 

p + (y" - t)p 

(pr — pr)xy"' + (y" — t) 2 (pr — pr) + (y" — t) 3 (pr — pr) 



x(p + (y" - t)p) 



V)3 



These expressions are substituted into the symmetry condition (2.2), and powers of y" are 
equated to yield an over-determined system of nonlinear ODEs for p and r. These are 
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easily solved with the aid of the contact condition (3.3). There are two sets of solutions. 
Either 

(p,r) = (c,ct), c 2 = a, (3.4) 



or 



(p,r) = (ct,c), c = -a. (3.5) 

Re-writing these solutions in terms of the original variables, we obtain eight inequivalent 
discrete symmetries that form a group isomorphic to Z4 x Z2. The group generators are 

Ti : (x,y,y') >-> (ix,-y,iy'), 

T2 : (x,y,y') h-> (y',xy' - y,x). 

The subgroup generated by Ti consists of four discrete point symmetries, which can be 
found without having to consider contact symmetries [7] . The four inequivalent non-point 
contact symmetries are obtained from the point symmetries by composition with T2, which 
is the prolonged Legendre Transform. 

Generally speaking, the larger the dimension of C, the easier it is to find the discrete 
symmetries. Nevertheless, the contact condition makes it possible to solve the governing 
equations, even when N = I. Suppose that C is one-dimensional and the ODE is written 
in canonical coordinates, so that the continuous symmetries are generated by 

X = dy. 

Then Lemma 1 gives X = bX, where 6/0, and therefore 

x = f(x,y'), y = by + g(x,y'). 
The contact condition yields 

y' = h(x,y'), 

where 

9x = fxh - by', g y/ = f y ,h. (3.6) 

Equations (3.6) are compatible if and only if 

fxhy' - fy'h x = b. (3.7) 

They can be integrated, once / and h are known, to determine g up to an arbitrary 
constant, which may be set at any convenient value to factor out equivalence under the 
one-parameter group generated by X. 

Consider the general third order ODE admitting the group generated by X: 

y"' = u;(x,y',y"). 

Substituting 

x = f(x,y'), y' = h(x,y') 
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into the symmetry condition, and then equating powers of y" , yields an over-determined 
coupled system of nonlinear PDEs. This system is precisely as intractible as the problem 
of using the symmetry condition alone to find all point symmetries of 

y" = u{x,y,i/). 

However, with the aid of the contact condition (3.7), the problem simplifies considerably. 
To illustrate this, consider the ODE 

y'" = y" 3 sin^y (3.8) 

whose only Lie contact symmetries are those generated by X = d y . The symmetry condi- 
tion gives the over-determined (but complicated) system 



fxh-xx fxxhx — h x sin 



fy'h X x + 2f x h X y ! fxxhy' 2f X yih x — 3h x hyi SUl 

fxhy'y' + 2f y 'h xy ' — fy'y'h x — 2f xy 'hy' = 3h x hy, sin (J^J > 
fy'hyiyi - fy'y'hy, + (f x hyi - fy'h x ) sin = hy, sin 



/ 
h 

This system can be greatly simplified by using (3.7) and its differential consequences, 
which reduces the first three equations to 

fxX = fxy' = h X = hy'yl = 0. 

Combining this result with (3.7) and the remaining symmetry condition gives 
/ = a(x + 2niry'), h = y', a £{-1,1}, n € Z. 

After solving (3.6) for g(x, y') and setting g(0, 0) = to factor out the continuous sym- 
metries, we obtain the following result. The inequivalent discrete contact symmetries of 
(3.8) form a countably-infinite group, which is generated by 

Ti : (x,y,y) ^ (-x,-y,y), 

T2 : (x, y, y') ^ (x + 27n/, y + iry' 2 ,y'). 
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4. Discrete uniform contact symmetries 

For ODEs of order n > 3, the Lie contact symmetries can generally be found systematically, 
and C is finite-dimensional. Second-order ODEs have an infinite-dimensional Lie algebra 
of contact symmetry generators, but they cannot all be found unless the general solution 
of the ODE is known. However, some contact symmetries may be found with the aid 
of a suitable ansatz for Q. For example, the restriction Q y > y > = enables all Lie point 
symmetries to be found systematically. 

Other restrictions on Q are possible. For example, the set of all uniform contact 
symmetries of a given ODE of order n > 2 is a finite-dimensional Lie group [8]. Uniform 
contact symmetries are of the form 

x = $(x,y'), y = ky + Q{x,y>), y' = ^(x,y'), k € M\{0}, (4.1) 

where the contact condition requires that 

ky' + G X = Qy, = $y,V. (4.2) 

Many of the most commonly-occurring contact symmetries are uniform, including all con- 
tact symmetries of the ODEs that were used as examples in the previous section. The 
generators of uniform Lie contact symmetries have characteristic functions of the form 

Q = ciy + (f>(x,y'), ci € R. 

The Lie algebra of these generators can be found systematically, for a given ODE of order 
n > 2, by equating powers of y in the symmetry condition. The constructions leading to 
Lemma 1 can be repeated, restricting attention to uniform contact symmetries, to obtain 
the following. 

Lemma 2. Every uniform contact symmetry T of an ordinary differential equation of 
order n > 2 induces an automorphism of the Lie algebra, C, of generators of one-parameter 
local Lie groups of uniform contact symmetries of the differential equation. For each such 
T, there exists a constant non-singular matrix (bfj such that 

Xi = b\X h (4.3) 

This automorphism preserves all structure constants. 

This lemma can be used to find all discrete uniform contact symmetries in exactly the 
same way as Lemma 1 is used. To illustrate this, consider the ODE 

which has a four-dimensional Lie algebra of uniform contact symmetry generators. The 
structure constants are simplest in the basis 

X 1 = 2y'd x + y' 2 d y , 

X 2 = ^xd x + ^yd y - ^y'd yl , 



X 3 = 2(y'y 3 d x + 3(yT 2 0y, 
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X 4 = ^xd x + -yd y + -y'd yl . 



The only non-zero structure constants are 



c 12 



c 21 



-34 



" c 43 



and so the Lie algebra is isomorphic to a(l) © a(l). After simplifying B as far as possible, 
using the relations (2.16) and the adjoint action of the continuous group, we obtain two 
possibilities. Either 



B 



or 



B 



a 














1 




























1 








a 














1 


p 














1 









a,/3€ {-1,1}, 



(4.5) 



a,0 €{-!,!}■ 



(4.6) 



The inequivalent discrete uniform contact symmetries are calculated in the same way as 
previously. Lemma 2 is used, together with (4.1) and the contact condition (4.2), to obtain 
the most general form possible for a uniform contact symmetry of (4.4). The symmetry 
condition for the ODE is then used to determine which of the possible solutions actually are 
symmetries. We find that the group of inequivalent discrete uniform symmetries of (4.4) 
is isomorphic to Z2 x Z2 x Z2, and is generated by 

Ti : (x,y,y') >-> {-x,y,-y'), 
T2 : (x,y,y') h-> (x, -y, -y'), 
T 3 : (Z, y, y') -> (xy' 2 , 2xy' - y, (y'y 1 ) . 



The discrete point symmetries generated by T± and T2 can also be found directly from 
the Lie algebra of point symmetries, which is Span(X2,X4). These inequivalent discrete 
symmetries are derived from the matrix B in (4.5), with a = (3; they map each of the o(l) 
Lie subalgebras to itself. However T^, which is derived from (4.6), interchanges these two 
subalgebras. Actually, (4.4) is merely one representative of a whole class of second-order 



ODEs that have inequivalent uniform discrete symmetries Tj 
of the form 



1,2,3; these ODEs are 



my 



y — (m + l)xy' ' 



m € 
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5. Contact symmetries of a nonlinear PDE 



The method described in section 2 generalizes to PDEs without difficulty; the correspond- 
ing algorithm for point symmetries will be discussed elsewhere [9]. The basic steps are 
shown here, using the potential hyperbolic heat equation 

utt + u t = — (5.1) 



as an example. The Lie algebra of contact symmetry generators is five-dimensional, with 
a basis 



x 1 = d t , x 2 = d x , x 3 = d u , 

X 5 = -xd x + ud u + u t d Ut + 2u x d Ux . 



X, 



e l d u 



These Lie contact symmetries are actually point symmetries, but there is still the pos- 
sibility that some discrete contact symmetries may be non-point symmetries, cf. (3.1). 
Using (2.16) and the adjoint action of the continuous group to simplify B, we obtain two 
possibilities. Either 

B = diag{l,6,a,/3,1}, a,/3e{-l,l}, b € M\{0}, 



or 



B = 



10 

6 

a 

/3 

2 








-1 



a,/3e{-l,l}, /)GK\{0}. 



(5.2) 



(5.3) 



B 



If B is of the form (5.2), the equations analogous to (2.13), (2.14) are 

~~Xii X ± x Xxu" 

X 2 t X 2 x X 2 u 

X 3 i X 3 x X 3 u 

X^t X4X X^u 

X$t X5X X§u 

whose general solution is 

x = bx + c 2 \u x \ 



t = t + Cl, 



1 











1 











1 








e - £ 





—x 


u 




1 




x 1 


2 , 


u 



au + (a — (3e C1 )u t + c 3 \u x \^ , Cj € 



The contact condition and the symmetry condition reduce these further, to 

(t,x,u) = (t,ax,au), a £ {—1,1}. 

The remaining inequivalent discrete contact symmetries are obtained from (5.3) in a similar 
way. To summarize the results: the inequivalent discrete contact symmetries of (5.1) form 
a group isomorphic to Z 2 x Z 2 , which is generated by 

Ti : (t,x,u,u t ,u x ) h-> (t,-x,-u,-u t ,u x ), 

T 2 : (t,x,u,ut,u x ) h-> (t + ln\u x \, u + u t , x + ut(u x ) _1 , -u t {u x )~ l , (ux)" 1 ) . 

The symmetry T 2 was known previously [15], but it was not known that T 2 and Til^ are 
the only non-point contact symmetries (up to equivalence). The method outlined in the 
current paper enables the user to completely classify all contact symmetries of a given 
differential equation with a known non-trivial Lie algebra. 
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